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ABSTRACT 



We present the results of a new model calculation of the 7-ray background 
produced by unresolved blazars, using the second EGRET catalogue and taking 
account of flaring. These results are compared to the preliminary 7-ray back- 
ground spectrum reported recently by the EGRET team. We find that blazars 
can account for the entire extragalactic 7-ray background observed by EGRET. 
In addition the EGRET spectrum shows some evidence of a signature for the 
blazar background, a concavity in the spectrum first pointed out in our earlier 
paper. 



Subject headings: galaxies: active - gamma-rays: theory - quasars: general 



1. Introduction 

An isotropic, thus presumably extragalactic, component of the cosmic 7-ray background 
above 30 MeV was first detected by the SAS-2 satellite QFichtel, Simpson, and Thompson 



1978| , Thompson and Fichtel, 1982| ) and has more recently been measured by the EGRET 
experiment on the Compton Gamma Ray Observatory ([Fichtel, 1995| ). The EGRET team 
has also discovered that the only significant extragalactic point sources which produce 7-rays 



in this energy range are the catagory of objects called blazars ( [Lhompson et al., 1995| , permer 



and Gehrels, 19951) . These objects are active galactic nuclei whose jets are nearly pointing 
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along our line of sight, resulting in enhanced luminosity and variability within a narrow cone 
around the jet axis, due to the highly relativistic motion of the source region within the jet. 
It can also be shown that physical mechanisms for producing diffuse radiation, such as the 
annihilation of antimatter or of supersymmetry dark matter particles, cannot produce an 
isotropic 7-ray flux above a few hundred MeV of the magnitude observed flStecker, 1989 ). 



Blazars would thus appear to be the only present candidates for producing the observed 
extragalactic 7-radiation in this energy range. 

With the new measurement of the extragalactic 7-ray radiation background (EGRB) by 
EGRET ( [Fichtel, 1995| ), and an enhanced blazar data set in the second EGRET catalogue 
( |l'hompson et a/., 1995| ), we have been led to reconsider an earlier model of 7-ray background 
production by blazars ( |Stecker, Salamon, and Malkan, 1993| , |Salamon and Stecker, 19TM| ), 



which we have now expanded to include the effects of both flux and spectral variability of 
blazars due to flaring. A particularly important observation by EGRET is the apparent 
hardening of spectra of blazars when they are in a flare state, which plays a key role in our 
analysis. Although this model is simplistic, it provides an excellent fit to both the shape and 
amplitude of the EGRB, while simultaneously accounting for both the number of blazars 
seen by EGRET, and their observed flux distribution. 



The Basic Model 



Our earlier model of 7-ray background production by blazars ( [Stecker, Salamon, and 
Malkan, 1993| , |Salamon and Stecker, 19M|) made the basic assumption that blazars seen in 



7-rays above 100 MeV are also seen in the radio (2.7 to 5.0 GHz) as flat spectrum radio 
quasars (FSRQs), with fewer blazars seen than FSRQs due to the existence of a narrower 
beaming angle in 7-rays than in the radio (see also Dermer, 1995). With the additional 
assumption that the radio and 7-ray luminosity of these objects are linearly related (see 
also Padovani, 1993), one is able to use pre-existing radio luminosity functions (RLF) to 
estimate the blazar 7-ray luminosity function (GLF), which can then be used to calculate 
the contribution of unresolved blazars to the isotropic 7-ray background. A constraint on the 
model is that the number of observed blazars in 7-rays estimated with this GLF correspond 
to the number actually observed by EGRET, given their instrumental flux sensitivity. 



We use here the same RLF as in our earlier work ( punlop and Peacock, 199HD , 
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moving density p r are, respectively, WHz -1 sr _1 and Mpc _3 -(unit interval of log 10 P) -1 . P r is 
the differential radio luminosity (per steradian) measured at 2.7 GHz; the RLF is constructed 
assuming pure luminosity evolution, with Q = 1.0 and H = 50 kms -1 Mpc -1 . 

The GLF is then given by 

p^(P 7f ,z) =rjp r (P jf /K,z), (2) 

where rj = (9 1 /9 r ) 2 is the reduction of the number of 7-ray blazars compared to FSRQs by 
virtue of a hypothesized smaller beaming angle # 7 , and k is the proportionality constant 
between the differential luminosities, P 7 / = nP r , which we evaluate at the fiducial energies 
Ef = 100 MeV and Uf = 2.7 GHz. The parameter k is estimated from the distribution of 
ratios of measured 7-ray and radio fluxes of the EGRET blazars, where radio fluxes at 2.7 
GHz were obtained from |Wall and Peacock, 1985|, |Kuhr et at., 1981], and |Ledden and U'DeTT 



1985 



With this GLF, and appropriate values for 77 and k, an estimate of the contribution 
to the EGRB from blazars was made ( [Stecker, Salamon, and Malkan, 1993| , [Salamon and| 



Stcckcr, 1994 ) and found to be a significant fraction of the measured EGRB. In particular, 
we noted that the spectrum of the EGRB due to blazars would be softer at lower energies, 
and harder at higher energies, a feature which appears to be present in the most recent 



EGRB measurement ( [b'ichtel, 1995 ). We also note that a recent analysis using EGRET data 



to obtain an estimate of the GLF at the high end of the blazar luminosity distribution is 
consistent with a luminosity and redshift dependence of source distribution given by Eq.[l] 
flChiang et ag 1995|) . 

An obvious oversimplification of the earlier model was its lack of accounting for the well- 
known time variability of blazars, a complex phenomenon which does not lead itself easily 
to simple characterizations. It is known that blazars can vary significantly over timescales 
of the order of a couple of days or less (|Jang and Miller, 1995| , [Wagner and Witzel, 1995| ). 



Significant flaring in the GeV 7-ray range is perhaps best illustrated by the intense flare 
which occurred during June 1991 in the source 3C279 ( Knitfcn et al, 19*93 ). Such flaring, in 7- 



rays as well as in other regions of the electromagnetic spectrum, appears to be associated with 
shocks propagating down relativistic jets and with the subsequent acceleration of relativistic 
particles by these shocks ( [Jang and Miller, 1995 , Valtaoja and Tcrasranta, 1995 ). 



In this paper we incorporate this phenomenon in a simple way, postulating the existence 
of two "states" of a blazar, a "quiescent" state and a "flaring" state, where the quiescent 
blazar luminosity is amplified by a fixed factor A during a flare. It is assumed that a blazar 
spends a fraction ( of its time in a flare state, where ( 1. 

In addition, the distribution of spectral indices of flare-state blazars is taken to be harder 
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Fig. 1. — A plot of the shift in measured blazar spectral index versus flux ratio for various 
blazars observed by ERGET during different viewing periods. The flux values used in the 
ratios were those obtained for the viewing periods in which spectral indices are determined. 
The datum for 0528+134 is taken from Mukherjee et a/.(1995). 



than that of quiescent-state blazars ( fSreekumar et al., 1995 ). This is supported by EGRET's 
observations of flares in which the hardening of the spectral index is statistically significant 
e.g., 0528+134 flMukherjee et al, 1995| ), and also as a statistical trend as seen in Figure 1. 
Here we have used the 2 nd EGRET catalogue ( Thompson et al, 1995|) and considered all 
blazars for which more than one spectral index is given (corresponding to different viewing 
periods). For each object, the ratio of the integral 7-ray fluxes of the two different periods 
is plotted against the spectral index shift, where we construct the flux ratio so that it is 
always greater than unity. A suggestive trend is apparent, indicating a typical hardening of 
the spectral index by ~ 0.3 during strong (A > 3) flares. 

The spectral index a is defined by Ej^(E) oc E~ a , where j^(E) is the differential 
number flux of 7-rays at energy E; equivalently, the 7-ray differential luminosity at the 
source varies with 7-ray energy as P^{E) = P 1 f(E/Ef)~ a , consistent with all observed 
blazar 7-ray spectra. The blazars' distribution of spectral indices a is determined from the 
2 nd EGRET catalogue, 

1 N 1 

#) = 4e^" MiW . (3) 
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where and cij are the catalogue's spectral index and error for the i th of N blazars. We 
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simply shift this distribution to obtain the separate distributions for the flare and quiescent 
states: pf(a) = p(a — Actf), and p q (a) = p(a — Aa q ), where we take Aa q = 0.20 and 
Aaj = —0.05, as discussed below. As there is no evidence for a change in the blazar spectral 
index distribution with redshift z in the EGRET catalogue, we assume none. 

Our GLFs for the two states of blazars, quiescent and flaring, is a modification of Eq.[| 

p 7 (P 7/ , Z) = (l- CVr(^, Z) + CVPrfa, Z), (4) 

K AK 

where the quiescent and flare 7-ray luminosities, P^ and P?, are related to the radio lumi- 
nosity as P!!j = kP t and P? = AnP r . We note that although variability is seen in the radio 
at higher frequencies, very little is observed at 2.7 GHz, the frequency at which our chosen 
RLF is determined ( [Reich, et al, 19931 , [Valtaoja and Teriisranta, 19951 ). Thus we need not 
be concerned with the effects of variability in our reference RLF. 



3. The Observed Flux Distribution 



For a blazar at redshift z with luminosity P y f (which is the luminosity evaluated at the 
fiducial energy Ef at the source) one can show that the integral flux F(E/) at Earth is given 
by 

F(E f ) = P, f /[a(l + zr +1 r 2 Rl], (5) 

where r is the co-moving coordinate of the blazar, and Rq is the present cosmological scale 
factor, R r = (2c/H ) 1 - (1 + zy 1/2 for fi = 1 and A = 0. (The absence of An in the 
denominator is because P 7 j is a differential luminosity per steradian.) 

The number of sources at redshift z seen at Earth with an integral flux F is given by 

dN 



dFdz 



Az AF = AnRy Ar p 7 (P 7/ , z) A(log 10 P 
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where the relation between F and P 7 / is fixed by Eq]5|. With this, the number of sources 
present per unit logarithm of flux F is found by integrating EqM over redshift, 



dN 



d(log 10 F) HI jo 
where we take z max = 5. 
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For the choice of parameters discussed in Sec.f| below, the number of blazars in the sky 
as a function of integrated flux at Earth is shown in Figure 2. Here we show the number 
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seen as quiescent-state blazars, and the number seen as flare-state blazars, and their total. 
Note that the two separate contributions represent the same underlying (radio) distribution, 
shifted along the abscissa by the flare factor A (A = 5), and along the ordinate by the flare 
duty cycle, C (C = 0.03). Also shown on the figure are the number of blazars seen by EGRET 
as a function of flux F(Ef = 100 MeV), binned logarithmically. We take as our sample the 50 
blazars identified in the 2nd EGRET catalogue, along with 14 unidentified EGRET sources 
above 30 degrees galactic latitude which are most likely to be blazars flMukherjee et aU, 
|1995|) . The agreement is seen to be quite good, both in the shape and amplitude of the 
experimental distribution. The dropoff of data below F = 2 x 10~ 7 cm _2 -s _1 is due to the 
loss of blazar detection efficiency near EGRET's minimum flux sensitivity of ~ 1 x 10~ 7 
cm~ 2 -s~ 1 . In addition, there is incompleteness in the sky survey at this sensitivity near the 
galactic plane, owing to the enhanced galactic foreground emission. 




observed integral flux (>100 MeV) (10 cm -s ) 



Fig. 2. — The the predicted distibution of sources at various flux levels using the model 
parameters described in the text. The open circles and their statistical error bars represent 
the 50 objects identified as blazars in the second EGRET catalogue, plus the 14 unidentified, 
high-latitude sources which are presumably blazars as well. The open triangles (no error bars) 
show only the 50 identified blazars' count per logarithmic bin. Note that the ordinate is the 
full sky count per one- fifth decade of flux. 



The curves in Figure 2 give the number of sources vs. flux predicted at a single point 
in time. At any one point in time, a certain fraction of the observed blazars will be in a 
flare state. At another point in time, however, a different set of blazars will be in a flare 
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state. Since EGRET has viewed the entire sky multiple times (although not uniformly after 
its inital full sky map), one expects that the number of flaring blazars seen by EGRET 
would be greater than that indicated by the flare curve of Figure 2. How much this should 
increase the contribution of the flare component to the total EGRET blazar count is not 
clear, as many of the flaring blazars are also detected in the quiescent state ([von Montigny 
et al, 1995|) . Complicating this issue is the fact that many of the multiple observations are 
made significantly off EGRET's vertical axis, increasing the minimum flux sensitivity for 
detection. The mean number of observations of an object (with and without detections) is 
~ 7, at various levels of sensitivity. Were we to increase the flare contribution to the total 
EGRET count by a factor of 3 to account for multiple viewings, we would still be reasonably 
consistent with the data (allowing for off-axis inefficiencies at the lower flux levels). 

An important point is that the distribution of spectral indices a measured by EGRET 
contains those from both the quiescent and flare state populations, with the flare states con- 
tributing heavily to the EGRET distribution, as Figure 2 and the above discussion indicates. 
Therefore the quiescent state population, which we show below dominates the extragalactic 
7-ray background radiation, has a softer distribution of spectral indices than that seen by 
EGRET. 



4. The Extragalactic Gamma-Ray background 

Unresolved blazars, i.e., those which are below detection threshold, will contribute to 
a diffuse, isotropic background of 7-ray radiation. We now consider whether unresolved 
blazars may be responsible for the entire observed extragalactic 7-ray background. 

The differential number flux ^ 1 of 7-rays received at Earth from a single source at 
redshift z is given by |Salamon and Stecker, 1994| 

f\ Eo , ,) = „ (8) 

dE Ac 2 E f {l + z) a+1 [1 - (1 + z)' 1 / 2 } 2 

where E is the 7-ray energy at Earth, Ef is the fiducial energy (100 MeV), and a is the 
spectral index of the single source. We note that we neglect the minor effect of 7-ray 
extinction off the IR-to-UV diffuse background ( [Stecker, deJager, and Salamon, 1992|) for 
our energy region of interest, 0.1 to 100 GeV. 

The total differential flux from unresolved blazars at a redshift z is found by integrating 
Eq|| over the GLF p 7 (P 7 /, z): 

dF f A/ma, dF^ 

— {E )=ATxRlr 2 drJdap{a)J p ^ — (E ,z)p y (R r f,z)dQog l0 F rf ), (9) 
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where -P 7 /, m ax is the luminosity of a blazar at redshift z which is just below the threshold for 
detection (found by subtituting EGRET's integral flux sensitivity F^ RET (100 MeV) into 
Eq.|5|), and an integration over the blazar spectral index distribution is made. 

We next integrate over redshift to obtain the total diffential flux at Earth at energy Eq. 
We substitute Eq.f| into Eq.|, expressing the luminosity integral as one over the RLF (Eq.|l]), 
obtaining 

§ (JSb) = h^To {/ So— (TT#W, J^T W (1 - QMPr, *) dPr 

+ fP f (*) IS*" (TT#W^ J]fc W C^H Z) dP r } . (10) 

The radio luminosity limits are given by P? >max {z) = P 7 /, max (z) / k, and P r f max (z) = P~,f, max {z)/AK. 

This equation gives the results shown in Figure 3, plotted as£ 2 ^ (the " energy flux" ) vs. 
E, for the parameter values k = 4 x 10~ n , r) = 1.0, £ = 0.03, A = 5, A« 9 = 0.20, and Act/ = 
—0.05. Both the amplitude and shape of the calculated energy flux E 2 ^ matches that of the 
data, with the exception of the two end data points, which have large systematic uncertainties 
(not shown) (P. Sreekumar, personal communication, 1995). We note in particular the role 
of the spectral index shift between the quiescent and flare state populations: if Eq.[l(] were to 
be integrated over a using just EGRET's observed spectral index distribution, the minimum 
in E 2 j^ would displaced to a much lower energy, in contradiction to the data. 

Thus we find that in order to match the observed spectral shape of the EGRB, we must 
invoke a hardening of the flare-state spectra with respect to the quiescent-state spectra, a 
feature which is already suggested by EGRET observation (Figure 1). In addition, we find 
that the entire EGRB can be explained as being due to the 7-ray emissions of unresolved 
blazars. 



5. Discussion 

We find that those unresolved blazars which make up the 7-ray background are primarily 
in the quiescent state, and therefore have a mixture of spectral indices shifted somewhat to 
the steep side (higher a) than those which make up the set of observed blazars. This is 
because the observed blazars contain a higher mixture of flaring sources than those which 
make up the background, owing to the high EGRET detector threshold and the shape of 
the high end of the GLF (Figure 2). Thus the EGRET blazar spectral index distribution 
will be intermediate between that of the flare and quiescent populations, with the mean 
quiescent- state index being Aa q = 0.20 higher than the observed EGRET distribution, and 
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Fig. 3. — The predicted background spectra from summing unresolved blazars compared 
with the preliminary EGRET data on the cosmic 7-ray background. The lower, middle, 
and upper curves show the contribution to the EGRB from flaring sources, from quiescent 
sources, and the their total. The curves shown are for the parameters indicated in the text, 
which are the same as those used to generate Figure 2. 



the mean flare-state index being |Aaf| = 0.05 lower than the EGRET distribution. The 
difference in mean spectral index between the flare and quiescent states, 0.25, is consistent 
with the trend indicated in Figure 1. In order to match the concavity of the observed EGRB 
spectrum, we require a spectral index distribution that is softer than that seen by EGRET, a 
condition that is met by this two-population model. We note that a general feature of models 
which posit that the EGRB is due solely to blazars is that the predicted EGRB spectra is 
concave and not a simple power law. This is due to the EGRB being a superposition of 
blazar power-law spectra of different spectral indices, where the softer components dominate 
at lower energies, and harder components dominate at the higher energies. 

The calculated EGRB spectrum in Figure 3 is shown down to an energy of 10 MeV. 
Of the few EGRET sources which have been detected by OSSE and COMPTEL below 
this energy, 7 by COMPTEL and 5 by OSSE ( McNaron-Brown et «/., 1995| and references 



therein), almost all show a flattening below ~ 10 MeV. There are not enough sources detected 
at the lower 7-ray energies to calculate a background, although it is reasonable to assume that 
the background from unresolved blazars should also flatten below 10 MeV. Note also that 
the calculated energy flux of the EGRB increases by more than a factor of 2 as one goes from 



-10- 



30 MeV down to 10 MeV energy. A preliminary analysis of COMPTEL data ( Kappadath ei 



U., 1995| ) gives a value for the GRB at 10 MeV which is more than a factor of 2 below our 



results. It is possible that the discrepancy here may be due to an oversubtraction of detector 
backgrounds in obtaining the residual GRB from the COMPTEL data, this analysis also 



having given much lower values than those obtained by APOLLO 15 flTrombka et ai, 1977|) 
at lower energies. 

Finally, we discuss the parameter values chosen to match the EGRET data. A single 
flaring factor of A = 5 is a gross simplification to what is obviously a continuous distribution 
of flaring factors, and is perhaps as large a value as one should take in this model, although 
a large A helps distinguish the separate roles of the two blazar states. A surprising result of 
this calculation is that we were driven to a large value of rj, which is the ratio of solid angles 
of the beamed 7-ray emission to the beamed radio emission; in fact, our optimal value is 
unity. This contrasts with the results of our earlier work ( ]5alamon and Stecker, 1994 ) and 



the model of Dermer (1995), where the emitted 7-rays are emitted into a much smaller 
solid angle around the jet axis compared to the beamed radio emission. This smaller 7-ray 
beaming angle was invoked in our earlier model to explain geometrically why fewer blazars 
are seen in 7-rays than in the radio. However, in the present model, the reason fewer 7-ray 
blazars are seen is because the ratio k between the 7-ray and radio luminosities is lower. 

The relationship between the parameters rj and k can be seen in Eq.^. To match the 
observed number of sources per logarithmic bin of flux, the quantity 7]k v is fixed, where V 
is the local power law of the radio luminosity distribution, p r (P r ) oc P~ T '. When we include 
the effect of flare amplification, k is lowered relative to our earlier model; thus 77 becomes 
larger than in our earlier model. Owing to the fact that we have introduced two states of 
blazar activity, we have now have the freedom to raise rj so as to match the amplitude of 
the EGRB, while simultaneously respecting the observed number of sources. An indication 
that this added degree of freedom is not totally ad hoc is seen by the fact that we require 
a difference in the mean spectral indices of our two blazar populations which is consistent 
with that suggested by EGRET observations. 

Our conclusion then is that both the amplitude and shape of the extragalactic 7-ray 
background spectrum as measured by EGRET can be explained as being due solely to the 
7-ray emissions of unresolved blazars, if one assumes a simple relationship between the 7-ray 
and radio luminosities of these objects. A critical feature of this simple model is that at a 
given time, a blazar is in one of two states: a quiescent state or a flare state. 



We thank P. Sreekumar for many helpful conversations and criticisms of the manuscript. 



- 11 - 



REFERENCES 

Becker, et al, 1991, Astrophys. J. Suppl. 75, 1. 

Chiang, J., C.E. Fichtel, C. von Montigny, P.L.Nolan and V. Petrosian, 1995, Astrophys. J. 
452, 156. 

Dermer, CD. and Gehrels, N., 1995, Astrophys. J. 447, 103. 
Dermer, CD., 1995, Astrophys. J. 446, L63. 

Dunlop, J.S. and Peacock, J.A. 1990, Mon. Not. R. Astron. Soc.247, 19. 

Fichtel, C.E., 1995, Proc 3rd Compton Observatory Symposium, Astron Astrophys. Suppl. , 
in press. 

Fichtel, C.E., Simpson, G.A. and Thompson, D.J. 1978, Astrophys. J. , 222, 833. 

Jang, M. and H.R. Miller 1995, Astrophys. J. 452, 582. 

Kappadath, et al. 1995, in Proc. 24th Int. Cosmic Ray Conf., Rome 2, 230. 

Kniffen, D.A., et al, 1993, Astrophys. J. All, 133. 

Kiihr, H., et al, 1981 , Astron. Astrophys. Suppl. 45, 367. 

Ledden and O'Dell, 1985, Astron. J. 298, 630. 

McNaron-Brown, K., et al. , 1995, Astrophys. J. 451, 575. 

Mukherjee, R., et al, 1995 Astrophys. J. 441, L61. 

Mukherjee, R., et al. , 1995, submitted to Astrophys. J. . 

Padovani, P. Ghisellini, G. Fabian, A.C and Celotti, A. 1993, Mon. Not. Royal Astron. Soc, 
260, L21. 

Reich, W. et al, 1993 Astron. Ap. 273, 65. 

Salamon, M.H. and F.W.Stecker, 1994, Astrophys. J. 430, L21. 

Sreekumar, P., et al. 1995, in preparation. 

Stecker, F.W. 1989, Nucl. Phys. B (proc. Suppl), 10B, 93. 

Stecker, F.W. De Jager, O.C and Salamon, M.H. 1992, Astrophys. J. 390, L49. 



Stecker, F.W., Salamon, M.H. and Malkan, 1993, Ap.J. Letters, 410, L71. 

Thompson, D.J. and Fichtel, C.E., 1982, Astron. Ap. 109, 352. 

Thompson, D.J., et al, 1995, Astrophys. J. Suppl, in press. 

Trombka, J.I., et al, 1977, Astrophys. J. 212, 925. 

Valtaoja, E. and Terasranta, H. 1995, Astron. and Ap. 297, L13. 

von Montigny, et al, 1995, Astrophys. J. 440, 525. 

Wagner, S.J. and A. Witzel 1995, Ann. Rev. Astron. Ap. 33, 163. 

Wall, J.V. and Peacock, J. A. 1985, MNRAS, 216, 173. 



This preprint was prepared with the AAS IATeX macros v4.0. 



